In[22]:=
m = IdentityMatrix[2]
l = 1
nl = {3.0132}
dl = {50*10^-9}
\[Lambda]l = {0.30184*10^-6}
\[CurlyEpsilon]l = {1/(4*\[Pi]*9*10^9*2.5752)}
For[i = 1, i < l + 1, i++,
 \[Delta] = 2*\[Pi]*nl[[i]]*dl[[i]]/\[Lambda]l[[i]];
 mi = {{Cos[\[Delta]], 
    I*nl[[i]]*(\[Mu]0*\[CurlyEpsilon]0)^0.5/\[CurlyEpsilon]l[[i]]*
     Sin[\[Delta]]},
   {I*\[CurlyEpsilon]l[[i]]/nl[[i]]/(\[Mu]0*\[CurlyEpsilon]0)^0.5*
     Sin[\[Delta]], Cos[\[Delta]]}};
 m = mi . m;
 print(m)]

PsTE = nl[[l]]*\[CurlyEpsilon]0/\[CurlyEpsilon]l[[l]]
PsTM = PsTE
P0TE = (\[CurlyEpsilon]0/\[Mu]0)^0.5
P0TM = 1
rTM = ((m[[1, 1]] + PsTM*m[[1, 2]])*
     P0TM - (m[[2, 1]] + PsTM*m[[2, 2]]))/
  ((m[[1, 1]] + PsTM*m[[1, 2]])*P0TM + (m[[2, 1]] + PsTM*m[[2, 2]]))
tTM = 2*P0TM/((m[[1, 1]] + PsTM*m[[1, 2]])*
      P0TM + (m[[2, 1]] + PsTM*m[[2, 2]]))
rTE = ((m[[1, 1]] + PsTE*m[[1, 2]])*
     P0TE - (m[[2, 1]] + PsTE*m[[2, 2]]))/
  ((m[[1, 1]] + PsTE*m[[1, 2]])*P0TE + (m[[2, 1]] + PsTE*m[[2, 2]]))
tTE = 2*P0TE/((m[[1, 1]] + PsTE*m[[1, 2]])*
      P0TE + (m[[2, 1]] + PsTE*m[[2, 2]]))
RTM = Re[rTM]^2
TTM = Re[tTM]^2
RTE = Re[rTE]^2
TTE = Re[tTE]^2


Out[22]= {{1,0},{0,1}}
Out[23]= 1
Out[24]= {3.0132}
Out[25]= {1/20000000}
Out[26]= {3.0184*10^-7}
Out[27]= {3.4335*10^-12}
Out[30]= 8.77589*10^11 \[CurlyEpsilon]0
Out[31]= 8.77589*10^11 \[CurlyEpsilon]0
Out[32]= (\[CurlyEpsilon]0/\[Mu]0)^0.5
Out[33]= 1
Out[34]= ((-0.999985+0. I)+(8.77576*10^11+0. I) \[CurlyEpsilon]0-(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5))/((-0.999985+0. I)-(8.77576*10^11+0. I) \[CurlyEpsilon]0+(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5))
Out[35]= 2/((-0.999985+0. I)-(8.77576*10^11+0. I) \[CurlyEpsilon]0+(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5))
Out[36]= (0. +(8.77576*10^11+0. I) \[CurlyEpsilon]0-(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+(\[CurlyEpsilon]0/\[Mu]0)^0.5 ((-0.999985+0. I)+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5)))/(0. -(8.77576*10^11+0. I) \[CurlyEpsilon]0+(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+(\[CurlyEpsilon]0/\[Mu]0)^0.5 ((-0.999985+0. I)+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5)))
Out[37]= (2 (\[CurlyEpsilon]0/\[Mu]0)^0.5)/(0. -(8.77576*10^11+0. I) \[CurlyEpsilon]0+(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+(\[CurlyEpsilon]0/\[Mu]0)^0.5 ((-0.999985+0. I)+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5)))
Out[38]= Re[((-0.999985+0. I)+(8.77576*10^11+0. I) \[CurlyEpsilon]0-(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5))/((-0.999985+0. I)-(8.77576*10^11+0. I) \[CurlyEpsilon]0+(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5))]^2
Out[39]= 4 Re[1/((-0.999985+0. I)-(8.77576*10^11+0. I) \[CurlyEpsilon]0+(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5))]^2
Out[40]= Re[(0. +(8.77576*10^11+0. I) \[CurlyEpsilon]0-(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+(\[CurlyEpsilon]0/\[Mu]0)^0.5 ((-0.999985+0. I)+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5)))/(0. -(8.77576*10^11+0. I) \[CurlyEpsilon]0+(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+(\[CurlyEpsilon]0/\[Mu]0)^0.5 ((-0.999985+0. I)+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5)))]^2
Out[41]= 4 Re[(\[CurlyEpsilon]0/\[Mu]0)^0.5/(0. -(8.77576*10^11+0. I) \[CurlyEpsilon]0+(0. +6.16713*10^-15 I)/(\[CurlyEpsilon]0 \[Mu]0)^0.5+(\[CurlyEpsilon]0/\[Mu]0)^0.5 ((-0.999985+0. I)+8.77589*10^11 \[CurlyEpsilon]0 (0. +(0. +4.74969*10^9 I) (\[CurlyEpsilon]0 \[Mu]0)^0.5)))]^2